[1] 



o 
o 

(N 

c 

a 

Hi 

in 



> 

o 
o 



Non-Gibbsianness of the invariant measures 
of non-reversible cellular automata with totally asymmetric noise 

Roberto Fernandez^ and Andre ToomQ 

Abstract 

We present a class of random cellular automata with multiple invariant measures which are all 
non-Gibbsian. The automata have configuration space {0, 1} z , with d > 1 , and they are noisy 
versions of automata with the "eroder property" . The noise is totally asymmetric in the sense that 
it allows random flippings of "0" into "1" but not the converse. We prove that all invariant measures 
assign to the event "a sphere with a large radius L is filled with ones" a probability [i l that is too 
large for the measure to be Gibbsian. For example, for the NEC automaton (— ln/^) x L while 
for any Gibbs measure the corresponding value is xL 2 . 
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1 Introduction 



Studies of cellular automata and of their continuous-time counterpart, the spin- flip 
dynamics, have been successful in determining how many invariant measures the au- 
tomaton or dynamics have. Much less is known about properties of these measures. A 
natural question is whether they are Gibbsian, that is whether they could correspond to 
measures describing the equilibrium state of some statistical mechanical system. There 
are two categories of evolutions — both with local and strictly positive updating rates— 
for which the answer is known to be positive: (1) If the updating prescription has a high 
level of stochasticity — high noise regime — , in which case Gibbsianness comes together 
with uniqueness of the invariant measure |XBj, [T^, [TBj ; and (2) if the updating satisfies 
a detailed balance condition for some Boltzmann-Gibbs weights [|20f[ . Known cases of 
non-Gibbsianness, on the other hand, refer to automata where the updating rates are 
either non-strictly positive [jT6j] , [[50], Chapter 7] or non-local Q23] . 



In this paper we present some examples of stochastic non-reversible automata — that is, 
automata not satisfying any form of detailed balance — , with multiple invariant mea- 
sures, all of them non-Gibbsian. Our class of automata can be seen as a generalization 
of the North-East-Center (NEC) majority model introduced in [££| and discussed in 
many papers. Its non-ergodicity was first proved in [|28| (see also the discussion in [|T5fl ) 



[3] 

and later by another method in [f2j. Also it was simulated more than once [[I], £2T], |22 . 
Models of this sort are obtained by superimposing stochastic errors (noise) to deter- 
ministic automata having the so-called eroder property: finite islands of aligned spins, 
within a sea of spins aligned in the opposite direction, disappear in a finite time. 

We allow only one-sided noise or stochastic error — a "0" can stochastically be turned 
into a "1" , but not the reverse. Thus some of our transition rates are zeros and therefore 
the "dichotomy" result of J2D], Corollary 1] is not applicable. Our work does not settle 



the long-standing issue of the Gibbsianness of the invariant measures of NEC models 
with non totally asymmetric noise. There are conflicting arguments and evidences 
for the model with symmetric noise: An interesting heuristic argument has been put 
forward [J3t], Chapter 5] pointing in the direction of Gibbsianness, and a couple of pioneer 
numerical studies yielded findings respectively consistent with Gibbsianness Q2Tj] and 



non-Gibbsianness |22jj . However, we hope that the simple non-Gibbsianness mechanism 



clearly illustrated by our examples could be a useful guide and reference for the study 
of the more involved two-way-noise situation. 

In our examples, non-Gibbsianness shows up in the same way as in the basic voter 
model Rlfifl : Large droplets of aligned ("unanimous") spins have too large probability 



for the invariant measures to be Gibbsian. More precisely, we show that once a suitable 
"spider" of "1" appears, the dynamics causes the alignment of the spins in a neighboring 
sphere. This sort of damage-spreading property (or error-correcting deficiency) implies 
that the presence of a sphere of "1" is penalized by the invariant measures only as a sub- 
volume exponential. This contradicts well known Gibbsian properties. In fact, we can 



[4] 



be more precise. Gibbsian measures are characterized by two properties uniform 



non-nullness and quasilocality. As we comment in Section |3|, the large probability of 
aligned droplets means that the invariant measures can not be uniformly non-null. More 
generally, such invariant measures can not be the result of block renormalizations of 
non-null, in particular Gibbsian, measures. Furthermore, known arguments |7j (briefly 
reviewed in Section |3] below), imply that if one of these measures is not a product 
measure, then its non-Gibbsianness is preserved by further single-site renormalization 
transformations. 



2 Simple examples 

Before plunging into the technical and notational details needed to describe our results in 
full generality, we would like to present some simple examples that contain the essential 
ideas. The examples are defined on the configuration space {0, 1} z . 



Example 1: The NEC model . Its deterministic version is defined by a translation- 
invariant parallel updating defined by the rule 

a&ftO, 0) = major{^(0, 1) , x\l, 0) , a;*(0, 0)} , (1) 

where denotes the configuration at site £ 2 2 immediately after the 

t-th iteration of the transformation and major : {0, l} 2fc+1 -h. {0, 1} is the majority 
function, i.e. the Boolean function of any odd number of arguments, which equals "1" 
if and only if most of its arguments equal "1". This prescription yields an evolution, 
which is symmetric with respect to the flip <->■ 1 [a function with this property is 
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called a self-spin-flip function in Section [| below]. We consider a noisy version, where 
in addition spins "0" flip into "1" independently with a certain probability e , while 
spins "1" remain unaltered. This corresponds to an stochastic updating 

Prob(x m (*, j) = 0\x t ) = (1 - e) [1 - x$£(i,j)] . (2) 

The "all-ones" delta-measure Si is invariant for this automaton. For small e there is 
at least another invariant measure (see lemma fO] below). 

Let us start with the following simple observations which are immediate consequences 
of the NEC rule ([I]) and the one-sidedness of the noise: 

(i) Horizontal lines (parallel to axis i ) filled with spins "1" remain invariant under 
the evolution. 

(ii) The same invariance holds for vertical lines (parallel to axis j ) filled with spins 

(iii) After one evolution-step (that is, after one parallel updating of all the spins), a line 
of slope —1 filled with spins "1" moves into the parallel line immediately to the 
South- West. 

(iv) If the (infinite) "spider" formed by the i-axis, the j -axis and the line i + j = is 
filled with "1", then after t steps the evolution causes the whole triangle '■ 
i,j < 0, i + j > -t] to be filled with "1". 

The last observation can be visualized as a displacement, at speed 1, of the "front" 
formed by the line i + j = , with a simultaneous displacement (here a trivial one), 
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at speed 0, of the "fronts" formed by the i - and j -axis. This combined displacement 
produces a growing triangle full of "1" . 

The same observations hold if full lines are replaced by finite segments, except that, 
depending on the values of neighboring spins, in each iteration each segment can lose 
one or both of the "1" at its endpoints. We conclude that if at some time the spider 

SP ( o,o),L = {(«, 0) G Z 2 : -8L < i < 4L} U 

{(0,j) G Z 2 : -8L < j < 4L}U{(», j) € Z 2 : i + j = , -6L < i < 6L} (3) 

is filled with "1", then after 4L iterations the "1" fill a triangular region that contains 
the sphere S(-l-l),l • Therefore, if /i is a invariant measure, 

Mls^^J > M1sp m ,) > £ 3(12i+1) • (4) 

We have denoted 1a , for Ac Z 2 , the event {x : x(i,j) = 1, (i, j) G A} . The last 
inequality in ((|) follows from the fact that a "1" has a probability at least e to appear 
at a given site because of the noise. As commented in Section |3], such a probability is 
too large for the invariant measure to be Gibbsian, or block-transformed Gibbsian. 

Example 2: North-South maximum of minima (NSMM) . The initial deterministic pre- 
scription is defined by 

4+1(0, 0) = = max{min(x t (0, 0) , x\l, 0)) , min(^(0, 1) , x\l, 1))} (5) 

plus translation-invariance. The corresponding evolution is not symmetric under flip- 
ping, unlike the previous example. The stochastic version is obtained by adding one- 
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sided noise as in (g). For small e this automaton has more than one invariant measure 
(see lemma |4.2j) . One of them is, of course, the "all-ones" delta-measure 5i . 



The mechanism for non-Gibbsianness for this model is even simpler to describe than for 
the NEC model. Indeed, it suffices to observe that whenever a horizontal line is filled 
with "1", then in the next iteration these "1" survive and in addition the parallel line 
immediately to the South becomes also filled with "1" . The same phenomenon happens 
for finite horizontal segments, except that each creation of a new segment filled with 
"1" can be accompanied by shrinkages of up to two sites (the spins at the endpoints) 
of all the previously created segments. We conclude that if the "spider" (which looks 
more like a snake in this case) 

SP(o,o),l = {(i, 0) G Z 2 : -3L < i < 3L} (6) 

is filled with "1" at some instant, then 2L instants later the "1" will cover at least a 
square region that includes the sphere SVo,-l),i ■ Arguing as for (£|), we obtain for all 
invariant measures fi the bound 

M (l V _ t) J > MlsiwJ > ^ L+1 > (?) 
which implies that fi is neither Gibbsian nor block-transformed Gibbsian. 

A comment by A. van Enter (private communication) gives a colorful description of the 
mechanism acting in both preceding examples: "the spider fills his stomach faster ( x L 
sites at a time) than his legs shrink ( x 1 sites at a time)". 
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Example 3: A non-example . The automata defined by the deterministic prescription 
4^(0,0) = 

major{min(y(0,2), x\-l,2)), min(x'(2, 0), x t {2, -1)), min(x*(0, -1), a;*(-l,0))} 

(8) 



followed by one-sided noise (g), also has multiple invariant measures; this follows from 
lemma ( |4.4 ) (see below) because its <7q is empty. Nevertheless, neither the mechanism 



of Example 1 (travelling fronts), nor that of Example 2 (growing strips) are present, so 
the theory of the present paper does not apply. 

3 Non-nullness and the probability of aligned spheres 

We present in this section the key property used in our paper to detect non-Gibbsianness. 
To state it in its natural generality we introduce some definitions. 

We consider a general space of the form Q = S z where S is some finite set, equipped 
with the usual product a -algebra. For A C Z d and z G Q we denote za the cylinder 

2a = {x G Q : Xi = Zi, i G A} . (9) 

Definition 3.1 A measure /i in Q is said to have the alignment-suppression prop- 
erty (ASP) if for every configuration z G £1 

-ln M (2 A )HA| (10) 

for every finite set Ac Z d . 
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Here and in the sequel / -< g or g y f , for / and g positive functions means that 
there exists a constant C > such that / > Cg . 

All Gibbs measures have the ASP property, but many non-Gibbsian measures too. We 
construct now a general class of measures with this property by considering renormal- 
ized measures having suitable non-nullness features. For this we consider an auxiliary 
configuration space Qq = . The single-site space S can be very general, not 
necessarily finite or even compact. We assume that there is a a -algebra on S and 
consider the usual product Borel a -algebra on Qq . A renormalization transformation 
from Qq to Q is a probability kernel T{ ■ \ ■ ) from Q to Q . More explicitly, for each 
u G f^o , T(-\u) is a probability measure in £1 , and for each measurable event A 
of Q , T(A\ • ) is a measurable function on Qq . In words, T(A\u) is the probability 
that, given a configuration u G , the "renormalized" configuration is in A . This 
represents a general stochastic transformation while deterministic transformations are 
the special cases obtained via delta- like prescriptions T( • \u>) . A block-renormalization 
transformation is a transformation, for which probabilities factorize in the following 
sense: to every i G Z d there corresponds a finite set B(i) C Z d , called block, with 
the following properties: 

(i) If two points are far enough from each other, the corresponding blocks are disjoint. 
That is, there is a positive do such that if the distance between k,£ G Z rf is 
greater than oIq , then B(k) fl B(£) = ( do = 1 for the renormalization transfor- 
mations used in statistical mechanics, while do > 1 for common cellular-automata 
transformations) . 
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(ii) If z'i, . . . , ik are sites in Z rf , and ai, . . . , are values in .4 , then 

fe „ 

T({ar il = ai, . . . ,x ife = a k } \ u) = H T l} {{x l3 = aj} \ u) B{lj) ) . (11) 

.7=1 

Our notation indicates that the functions T i] {{xi j = aj} \ • ) depend only on the values 
of U£ for £ G B(ij) (i.e., they are measurable with respect to the a -algebra generated 
by the cylinders with base in B(ij) ). Examples of such transformations include decima- 
tion, (deterministic), Kadanoff transformations (stochastic), majority rule, sign fields 
and transitions of cellular automata (the last three can be deterministic or stochastic, 
depending on the setting). 

The kernel T naturally induces a transformation at the level of measures: Each proba- 
bility measure p on Qq is mapped into a probability measure pT on Q — the renor- 
malized measure — defined by 

j n f(x)( P T)(dx) = j n [j n f{x)T{dx\u)] P {du), (12) 

for all suitable / (e.g. continuous or non- negative measurable). For each mea- 
sure p on Q an d each block B(i) let us consider the conditional probabilities 
p(du>B(i) | oj x d \B(i)) ■ F° r a given transformation T we single out the set Vt of measures 
on Q,o that admit conditional probabilities such that 

if I* <4^ f BW / = ^ ' UB{ ^ ' p ( duJB ^ I u i d \B{i)) > 6 > ( 13 ) 

for some 5 > . We denote V the union of these families Vt over all block- 
renormalization transformations T . Here is our key characterization. 



Theorem 3.1 Every measure in V has the alignment-suppression property. 
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Proof. Let T , p be such that \i = pT . By property (ii) above, there exists a 
constant 7 > (proportional to do ) such that for any Ac 2 d there is a family of 
sites ii : . . . £ A with > 7|A| , all of which are far enough from each other and 
therefore the blocks B(ii), . . . , B(ik) are disjoint. We therefore have that for every 
z e Q 

p(T({x n = z n } I •) I Vtfwh)) II ^({^ = I UiHii)) P{du) 
< (1 " <*) / A V; ; ({.r = I u; sfe) ) p(do>) . (14) 

j=2 

This inequality is an immediate consequence of condition (|T3|) . After k iterations of 
this procedure we obtain 

m(*a) < (l-^) fc < (1-^) 7|A| . q.e.d. (15) 



The class V of measures is a very large class. It contains practically all block trans- 
formations of Gibbs measures with finite alphabet obtained via standard statistical 
mechanics prescriptions (decimation, Kadanoff, majority rule, etc), plus the measures 
generated by finite-time evolutions of usual cellular automata prescriptions. There is 
by now a vast literature about such measures — see, for instance, [j5], |TB|, |3f; for re- 
cent reviews with many references see H, [TO], [IT], [8j]- - showing that many of them are 
non-Gibbsian. In fact, the family Vi , where I is the identity, includes all uniformly 
non-null measures. These are measures fi that have, for each finite region Ac Z d , 
uniformly bounded conditional probabilities niduo^ \ cu 2 d \A.) > that * s ' sucn that there 
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exist 5a > with 

min inf h({xa = oa} | w z J > 5a- (16) 

We have denoted o-a = (&«)ieA • Gibbs measures are uniformly non-null — and in addi- 
tion quasilocal (the finite-volume conditional probabilities are continuous functions of 
the external conditions cJ z d A ) — hence they also belong to Vi . Property QI3D seems 
to be more general than usual non-nullness, in particular it does not depend on the 
existence of a whole system of conditional probabilities. 

The invariant measures of the automata of the present paper, on the other hand, do 
not have the alingment-suppression property, hence they do not belong to the class V . 
They therefore can be neither Gibbsian nor uniformly non-null nor block-transformed 
Gibbsian. As further examples of measures without the ASP we mention the invariant 
measures of the basic voter model [M , the invariant measure of some non-local dynamics 



23fl , an d the sign-fields of massless Gaussians |T4|, |5j, anharmonic crystals [|6], Section 



4.4] and SOS models g, [T7j]. 

For measures fi having a well defined relative entropy density s( • , the alignment- 
supression property ( |TD| ) implies that s(5 z \fi) > for every periodic configuration 
z G £1 . The relative entropy density is known to exist for translation-invariant Gibbs 
measures |T2|, Chapter 15]. Recent work in [jZ5| shows that it is also well defined for 



most translation-invariant measures obtained through block transformations of Gibbs 
measures. Because of this, the non-Gibbsianness resulting from the lack of ASP has 
often been interpreted as "large deviations probabilities that are too large" for Gibb- 
sianness. The non-Gibbsianness (non-nullness) criterion obtained by falsifying Theorem 
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however, is a more general argument that needs neither translation invariance of fi 
nor the existence of the entropy density. 

For completeness, we mention a further result obtained in [[7j. 

Theorem 3.2 Suppose \x is a measure in Q such that (i) it violates the ASP property 
for some periodic configuration zGO, and (ii) it is not a product measure. Then, for 
every single-site block-renormalization transformation T (i.e. a transformation defined 
by blocks B(i) formed by only one site), the measure fiT is not Gibbsian. 



This result follows from the fact that such a violation implies that s(5 z \fi) = , which 
in turns implies that s(5 z T\fiT) = . If fiT were Gibbs, by a well known result [|I2| , 



Theorem 15.37] the measure 5 Z T would be Gibbs for an equivalent interaction. But 
this impossible because the latter is a product measure and the former is not. Note 
that if T corresponds to a not-totally asymmetric noise, the measure fiT is uniformly 
non-null. Hence its non-Gibbsianness would correspond to lack of quasilocality. 

For the automata of this paper, we suspect that many of its invariant measures are 
non-product. 

4 General Results 



We now describe a large family of automata exhibiting a general version of the non- 
Gibbsianness mechanism of the first two examples in Section @. Throughout the 
article we consider the d -dimensional integer space 7L d with d > 1 embedded into 
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the d -dimensional real space IR d with the same axes and Euclidean norm II • II . The 
configuration space is Q, = {0, 1} z . We first need some definitions. 

For any i G Z d we denote Tj : Q — ► Q the translation of £1 defined by (r^ x)j = . 
Any function / : Q, — ► {0, 1} will be called a transition function. Given any transition 
function / , we define the corresponding operator Df : Q — > Q by the rule 

V*G Z d : (Dfx) i = f(T i x). (17) 
We call / : Q — > {0, 1} standard if it has the following three properties: 

1) / is local, i.e. there is a finite set A C Z d —the support of / - - such that 
f{x) = /(xa) • Given A , we denote p the maximum of ||i|| for i G A . 

2) / is monotonia, that is (V % : < =4> /(x) < f(y). 

3) / is not a constant. (Otherwise our theorem is either trivially true if / = 1 or 
trivially false if / = .) 

Since / is monotonic and non-constant, 

/("all zeros") = and /("all ones") = 1. (18) 

Let M. denote the set of probability measures on Q (on the a -algebra generated by 
cylinder sets). For any e G [0, 1] we define one-sided noise iV e : M. — > M. as follows: 
when applied to a measure 5 X concentrated in a configuration x = (xi) , it produces a 
product measure N £ 5 X , in which the % -th component equals 1 with a probability 1 if 
Xi = 1 and with a probability e if Xi = . 
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For any x G VL we denote its indicator Ind(x) = {i G 2 d | Xi = 1} . Conversely, for 
any 5 c Z d we denote Conf(S) that configuration, whose indicator is S . 

Let us call an element of JR d a direction if its norm equals 1. For any direction p we 
call a front with this direction any configuration whose indicator has the form 

{i G Z d | (i,p) < C}, (19) 

where C is a real number and (•, •} denotes the scalar product in JR d . It is evident 
that for any standard / the operator Df transforms any front (|K)|) into a front with 
the same direction, C being substituted by C + V P , where V p does not depend on C . 
We call V p the velocity of Df in the direction p . 

Let us call a configuration x G f2 invariant for if Df x = x . Given x, |/ 6 fl , we 
call y a finite deviation of x if the set of those i G Z d for which ?/j 7^ is finite. 
We say that an invariant configuration x attracts Df if for any its finite deviation y 
there is a time t such that Djy = x. 

Theorem 4.1 Take any standard f , such that "all ones" attracts Df , and make any 
one of the following two assumptions: 

a) V p + V-p > for all directions p . 

b) There is a direction p such that V p + V- p > . 

Then for any e > all the invariant measures of N £ Df satisfy 

-lnMlSoJ^- 1 . (20) 



If e = , our theorem may be false, for example if D is the identity. Notice also that in 
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the case b) our assumption that "all ones" attracts Dj is redundant because it follows 
from b). 

Let us present some further considerations that clarify the statement of the theorem. 
Given any non-constant amne function : lR d — > IR and two numbers C\ < C2 , we call 
a layer any configuration Conf{z G Z d C\ < (f)(i) < C2} ■ We call the thickness of this 
layer the distance between the hyperplanes (j) = C\ and (j) = C2 , that is (C2 — C\)/\4>\ , 
where | • | is the norm. We call a layer thick-enough if its thickness is not less than 2p . 

We call the two normal unit vectors to hyperplanes (j) = const the directions of this 
layer. If / is standard, Df transforms any thick-enough layer into a layer with the 
same directions, the thickness of the layer changing by V p + V~ p . The condition a) of 
our theorem means that thickness of any thick-enough layer does not decrease and the 
condition b) means that thickness of some layer increases under the action of Df . 

Of the examples of Section 0, the NEC automaton satisfies condition a), while the 
NSMM automaton satisfies condition b) for p = (0, 1) . For the non-example, however, 
V P + V-p < for all directions p . In all the three cases / [given, respectively, by (|Tj), 
and flSj)] is standard, and both "all zeros" and "all ones" attract Df . 



The NEC example is representative of a class of models with a further duality property. 
For any X{ G {0, 1} we denote ->Xi = 1 — X4 . Accordingly, if a: is a configuration, ->x 
is another configuration such that {~ , x)i = ->(xi) . Any transition function / has an 
associated spin-flip function denoted -1/ and defined by the identity ~>f(x) = ■ 
f\ Let us call / self-spin-flip if it coincides with its spin-flip. If / is standard and 

3 In the theory of Boolean functions ->/ is called dual, but in the theory of random processes the word "duality" is used 
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self-spin-flip, then V p + V- p = , so the thickness of all layers does not change under 
the action of Df . For example, the function major(-) , described above, is self-spin-flip. 



It is evident that under the hypothesis of Theorem |4.1|, the measure S\ is invariant for 



any superposition N £ Df . Hence, the theorem is not trivial only if the automata have 
more than one invariant measure. This is ensured by the following lemma. 

Given / , let us call a set Sc. 7L d a one-set if /(Conf(S)) = 1 . Since one-sets belong 
to Z d , they belong to IR rf , where we can consider their convex hulls, the intersection 
of which is denoted o\ . In the analogous way we call a set S C Z d a zero-set if 
/(Conf( Z d — S)) = and denote ctq the intersection of their convex hulls. 

Theorem 4.2 For any operator Df defined by (\T1[), where f is standard, the following 
four statements are equivalent: 

1) N £ Df has more than one invariant measure for some positive e . 

2) The configuration "all zeros" attracts Df . 

3) ctq is empty. 

4 ) There are a natural number m < d + 1 and m affine functions 

(pi,..., (j) m : IR d — > IR such that: 

i) for every j G [l,m] the set {p G Z d : (f)j(p) < 0} is a zero-set. 

ii) 0i + • • • + (fi m = const > . 

Hi) There is a rational point p G IR d such that 4>j{p) > for all j G [l,m] . 



for another purpose. 
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5 Proof of theorem |z£72 



If we omit the condition iii) in 4), our theorem |4.^| almost follows from theorems 5 and 6 



and lemma 12 of |j29(| . However, there is some difference, so for the reader's convenience 



we completely deduce 4) from 3). 

Suppose that gq is empty. Every zero-set can be represented as an intersection of 
several zero-half-spaces, i.e. half-spaces, which are zero-sets, where a half-space is a 
subset of JR d , where some non-constant affine function does not exceed zero. Thus 
there are several zero-half-spaces, whose intersection is empty. Everyone of them can 
be represented as {p G IR d | fi(p) < 0} , where are affine functions on IR d . We 
can choose these functions so that they have no common direction of recession (that 
is, no direction p such that fi(p) < fi{0) for all i), which allows us to apply to 
them Theorem 21.3 on page 189 of [129]. Since the intersection of our zero-half-spaces 



is empty, the case (a) of this theorem is excluded in the present situation, whence the 
case (b) takes place, which amounts to our conditions i) and ii) in 4), the products \fi 
mentioned in the case (b) serving as our fa . We may assume that our m is the 
minimal for which there are functions satisfying i) and ii). Based on this, let us prove 
statement iii) using the following lemma, which is a direct consequence of Theorem 21.1 
on page 186 of jZBj: 



Lemma 5.1 Let fa, . . . , <f) m be affine functions on R d . Then one and only one of the 
following alternatives holds: 
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(a) There exists some x G IR d such that 4>\{x) > 0, . . . , 4> m (x) > 0; 

(b) There exist non-negative real numbers X\, . . . , X m , not all zero, such that the sum 
Ai^i(x) + • • • + \ m (j) m (x) is a non-positive constant. 

Let us assume that the case (b) takes place in our situation. We may assume that A m 
is the greatest of Ai, . . . , A m , and therefore positive. From the statement ii) of 4), not 
all Xi are equal to X m . Let us divide all terms by A m : 

H h + 0m = const < 

and subtract this from the statement ii) of 4): 

f 1 - ^-W + • • • + ( 1 - ^] 0i = const > 0. 
V X m J \ X m ) 

Here all coemcients are non-negative and not all are zero. Therefore the functions 

(1 — Xi/X m )(f)i for i = 1, . . . , m — 1 also satisfy the conditions i) and ii) of 4) with a 

smaller value of m , which contradicts our assumption. Thus case (b) is excluded, so 

case (a ) takes place, whence there is a point p G IR d where all 4>j{p) > . Since all 

are continuous, there is a rational point with this property also, whence condition iii) 

of 4) follows. q.e.d. 

6 Proof of Theorem gTl] 

6.1 Proof of ([2D]) in case a) of the theorem 



Rewording Lemma |4.2| for the case when and 1 are permuted, we see that whenever 
/ is standard and "all ones" attracts Df , there exist a natural number m < d+ 1 and 



[20] 



m affine functions 0i,...,0 m : IR — ► R such that: 

•d . ± i 



i) for every j G [l,m] the set {i G Z' 

ii) 0i + • • • + (f> m = const > . 

,d 



i) < 0} is a one-set. 
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iii) There is a rational point p G H such that > for all 

j G [l,m] . 



For instance, for the NEC example there are m = 3 such affine functions, whose level 
lines are horizontal, vertical and lines of slope —1 respectively. 

For every j let us denote (f)j = — 4>j(0) whence 0j = <pj + 4>j(0) , where <pj is the 
linear part. Notice that |0j(O)| < \4>j\ ■ p and that 0i(O) + - • - + m (O) > . Notice also 
that if / is standard, "all ones" attracts Df and V p -\-V- p > for all directions p, then 
for any j G [l,m] and any thick-enough layer y = Confjz G Z rf I C\ < 4>j(i) < C2} 



lnd(Df y)^{i£ 7L d I d + 0,(0) < ^-(i) < C 2 + ^(0)}. 
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(See immediately above flI9|) for the definition of Conf and Ind.) 



Lemma 6.1 Take any standard f and assume that "all ones" attracts Df and that 
V P + V-p > for all directions p . Take x* defined by 



Ind(z*)= U Z d I |^(i)| < 2p-\cf> j \}. 

l<j<m 



(23) 



T/ien for t = 0, 1, 2, 3, . . . £/ie indicator of D l x* includes the union A t U B t , where 



A t = U )«e I !£;(«)-*■ 0j(O)| <2p- |^|} 



(24) 



1<J<TO 



and 



l<j<m 



(25) 
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[For the NEC example of Section [2j, this lemma corresponds to observation (iv).] 



Let us prove this lemma by induction. Base of induction: Since Aq coincides with 
Ind(x*) and Bo C Ao , our statement is true for t = . 

Induction step. Let us suppose that Ind(D t x*) D A t U B t , take any z G A t+ \ U 
and prove that i G Ind(Z) m x*) . Let us consider two cases. 

Case 1. Let i belong to A t+ \ . Then our statement follows from (|Z2|). 

Case 2. Let i belong to B t +\ , but not to A t+ i . Then 

^(i)-(t + l)-^(0)<-2p.|^| 

for all j G [l,m] . Notice that 

~4>i(i + v k ) < (f)j(i) + \(f)j\ ■ \v k \ < (f)j(i) + \(f)j\ ■ p. 

Therefore 

^{i + v^-t-^O) <0 j .(z) + |0 j |. /O -(t + i)0 j (o) + ^-(O) < -2p.|0il + l0il'P+0i(O) < 0. 
Thus 

j + AcB(C Ind(D* x*). 
Hence from (|T^) i G Ind(D i+1 x*) . Lemma |b."l| is proved. 



Lemma 6.2 Under the hypotheses of Lemma \b\ 1\, there is a positive constant a > 
such that for all t = 0, 1, 2, . . . the set B t defined by contains a sphere in ~ZL d 
with the radius a ■ t . 
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Proof. In fact we shall prove that 

V % G Z d , t = 0, 1, 2, . . . : \i + t ■ p\ < a -t => i G B t , 

where p is that rational point where all 4>j{p) > , whose existence is provided by hi). 
Let us denote k,j = 4>j{p) > and a = mirij , that is the minimal distance 
from p to the hyperplanes </> j = . Let us consider three cases. 

Case 1: p = . Then (ftj(0) = Kj > for all j . Now let us take any point i in the 
sphere with the radius a ■ t and center . This means that 

\i\ < a ■ t = mm(0j(O)/|0j|) • t. 

Then 

< \i\ ■ < 0j(O)/|0i| • * • = t ■ (f)j(0) 

for all j , whence i G B t . 

Case 2: p G 7L d . Then along with our operator Dj we consider another operator 
D g , where g{x) = f(r p x) . The function g is also standard, D g is also attracted by 
"all ones" and the affine functions provided for D g by iii) of (|2T|) can be obtained from 



those for Df by the same translation, so their values at are Ki, . . . , K m > , whence 
D g fits our case 1. So the set B t for D g contains a sphere with the center and radius 
a-t. Since Df commutates with all translations, the set B t for D g results from the 
set B t for Df by a translation at t-p. Thus the set B t for Df results from B t for 
D g by the opposite translation, whence it contains a sphere with the center —t ■ p and 
the same radius. 
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Case 3: p is any rational point. Let us denote q the least common denominator of all 
the coordinates of p and immerse our Z d into the set Z d , where Z q = {n/q n G 
Z } . Let us denote Q q = {0, 1} z « . Now / can be considered as a function g from 
Q q to {0,1} . Now let us "stretch" Z d q to turn it into Z d . Under this stretch the 
function g remains standard and "all ones" still attracts D g . In addition to that, the 
affine functions for D g with the properties Q2TQ now can be obtained from those for Df 
by a homothety with coefficient q . Therefore their values at the integer point q • p are 
Ki, . . . , K, m > . So D g fits our case 2, whence the set B t for D g contains a sphere 
with the center —t-q-p and radius a • q • t , whence the set B t for contains a 
sphere with the center —t ■ p and radius a • t . Lemma pT^ is proved. q.e.d. 



Now let us prove (\20j). From monotonicity it is sufficient to prove this inequality for 



H = (N £ DfY So for some t . Let us choose t\ such that a • t\ > R + d . Then, 
taking x* defined by (|23|) as the initial configuration, after t\ time-steps we obtain a 
configuration, whose indicator contains a sphere with the radius R + d and therefore 
contains a sphere with the radius R and center at some integer point p . However, 
what we actually need is a finite deviation from "all zeros", which coincides with x* 
only within a sphere with the radius R-\-t\ • p and has zeros outside it. The cardinality 
of its indicator does not exceed C{R d ^ 1 + l) with an appropriate C. Translating this 
configuration at the vector —p , we obtain another configuration, which fills with ones 
a sphere with radius R and center at the origin after t\ time-steps. The probability 
that the actual configuration's indicator contains this configuration is not less that 
£ c{R d ^ w h ence (jflj^ follows. q.e.d. 
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6.2 Proof of ([2D]) in case b) of the theorem 

This time we define x* as follows: 

Ind(x*) = {ie 7L d \ \(i,p)\ < p}. 

Then for all £ = 0,1,2,... 

Ind(L>}0 D {i G Z d | -p + t- V- p < (i,p} \ <p + t-V p } . 

Here the right side is a layer with the thickness 2p + t(V p + V- p ) . Given any R > , let 
us choose the minimal integer t\ for which 2p + t\(V p + > R+ d . Then indicator 
of D l j x* contains a sphere with an integer center and radius R . If we take an initial 
condition which coincides with x* within a sphere with the center at the origin and 
radius R+d+t\ -p , we shall obtain the same result. This configuration has C(R d ~ 1 -\-l) 
components that equal 1, where C is an appropriate constant. Further we argue like 
in case a). q.e.d. 

7 Final notes 

Note 1. Using minoration arguments, is is easy to expand our theorem to some random 
cellular automata, which cannot be represenred as N £ Dj . Using the same A as before 
and choosing transition probabilities 6{x\y&) for all x G {0,1} and y G {0, 1} A , we 
can define a random cellular automaton as an operator P : M. — > A4 which transforms 
any 5 y , where y G f2 , into a product-measure in which the probability that the i -th 
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component equals x is 0(x\yi + &) . This operator majorates N £ Df if 

f = 1 if f(n,A = 1 



= 1 if/(y A ) = l, 
>£ if /(i/a) = 0. 



As soon as this condition holds and Df satisfies conditions of our theorem, all invariant 
measures of P also satisfy Q2UQ and therefore are non-Gibbs. 



Note 2. In some cases it is possible to obtain a stronger estimation than (|2U|). Let 
d > a > and f{x) equal 

min max ■ ■ ■ , id) 

ii,...,i a e{0,l} z a+ i,...,z d e{0,l} 

where zi, . . . , are the coordinates of Z d . In this case 

-ln/i(l(5 0) L))^^ a , 

where ji is any invariant measure of N £ Df . If a < d — 1 , this estimation is stronger 
than ( |201 ) . This estimation can be proved in the same manner as in the case b), only 
x* now is defined by the condition: 

x* = 1 if max(|vfi|, . . . , < const. 

Note 3. Given a standard / , let us assume that "all zeros" attracts Df . Then we 
hope to estimate — ln//(l(So,z,)) from below as follows: 

-lnfi{l{S 0jL ))y L. 

If we succeed, this will settle the question of asymptotics of — ln//(l(5'o,i)) in some 
cases, e.g. in our examples 1 and 2, 
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Note 4. Those conditions under which our theorem holds and is non-trivial can be 
satisfied only for d > 1 . However, a statement similar to our theorem for the one- 
dimensional case was proved in j27j . Namely, it was proved that all non-trivial invariant 



measures of a class of one-dimensional random cellular automata did not belong to a 
class, which included all Markov measures. 
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